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Abstract 

Glueball spectrum is studied using an improved gluonic action on asymmetric lat- 
tices in the pure SU (3) gauge theory. The smallest spatial lattice spacing is about 
0.08/m which makes the extrapolation to the continuum limit more reliable. In 
particular, attention is paid to the scalar glueball mass which is known to have 
problems in the extrapolation. Converting our lattice results to physical units us- 
ing the scale set by the static quark potential, we obtain the following results for 
the glueball masses: Mg{Q'^~^) = 1730(90)Mey for the scalar glueball mass and 
Mg(2++) = 2400(95)Mey for the tensor glueball. 
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1 Introduction 



It is believed that QCD is the theory which describes strong interactions 
among quarks and gluons. A direct consequence of this is the existence of 
excitations of ghionic degrees of freedom, i.e. glueballs. However, due to their 
non-pcrturbativc nature, the spectrum of glueballs can only be investigated 
reliably with non-purterbative methods like lattice QCD [1-4]. Recently, it 
has become clear that such a calculation can be performed on a relatively 
coarse lattice using an improved gluonic action on asymmetric lattices [5-7]. 
In this paper, we present our results for a glueball spectrum calculation. The 
spatial lattice spacing in our simulations ranges from 0.08/m to 0.25 fm which 
enables us to extrapolate more reliably to the continuum limit. The improved 
gluonic action we used is the tadpole improved gluonic action on asymmetric 
lattices as described in [6,7]. It is given by: 

bTrPij 1 TrRij 1 TrRji 
9 Cuj ~ 36 ~ 36~^ 

In the above expression, (3 is related to the bare gauge coupling, ^ = ag/at is 
the (bare) aspect ratio of the asymmetric lattice with ag and a* being the lattice 
spacing in spatial and temporal direction respectively. The parameter Us is the 
tadpole improvement parameter to be determined self-consistently from the 
spatial plaquettes in the simulation. Pij and Pq* are the spatial and temporal 
plaquette variables. Rij designates the 2x1 Wilson loop (2 in direction i 
and 1 in direction j). Using spatially coarse and temporally fine lattices helps 
to enhance signals in the glueball correlation functions. Therefore, the bare 
aspect ratio is taken to be some value larger than one. In our simulation, 
we have used ,^ = 3 for our glueball calculation. It turns out that, using 
the non-pcrturbativcly determined tadpole improvement [8] parameter Ug, the 
renormalization effects of the aspect ratio is small [6,7], typically of the order 
of a few percent for practical values of P in the simulation. This could also 
be verified by measuring corresponding Wilson loops, which will directly yield 
the renormalized aspect ratio. For the moment, we will ignore their difference 
and simply use the bare value of It is also noticed that, without tadpole 
improvement, this renormalization effect could be as large as 30 percent [6,7]. 

This paper is organized as follows. In the next section, we will explain our 
calculation of the Wilson loops, glueball correlation functions which give us 
the mass values of the glueballs in various symmetry sectors. Finite volume 
errors are discussed and more importantly, finite lattice spacing errors are 
analyzed. Special attention is paid to the scalar sector where the extrapolation 
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used to be troublesome at coarse lattices. In the last section, we will have some 
discussion of out result and conclude. 



2 Monte Carlo Simulations 

Glueball spectrum calculations in pure gauge theory basically involve three 
stages. The first stage, gauge field configurations are generated using some 
algorithm. We have utilized a Hybrid Monte Carlo algorithm to update gauge 
field configurations. Several lattice sizes have been simulated and the detailed 
information can be found in Table. 1. For each lattice with fixed bare param- 

Table 1 

Simulation parameters and the corresponding lattice spacing in physical units ob- 
tained from Wilson loop measurements. 
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eters, order of 10^ configurations have been accumulated. Each gauge field 
configuration is separated from the previous one by several Hybrid Monte 
Carlo trajectories, typically 5 ~ 10, to make sure that they are sufficiently 
decorrelated. Further binning of the data has been performed and no notice- 
able remaining autocorrelation has been observed. 

The second stage of the calculation is to perform measurements of physical 
observables using the accumulated gauge field configurations. In fact, two in- 
dependent measurements have to be done. One is to set the scale in physical 
units, i.e. to determine the lattice spacing in physical units. This is neces- 
sary since there is no scale in a pure gauge theory. The second is to measure 
glueball mass values in lattice units. This is done by measuring various glue- 
ball correlation functions. With the scale being set in the first step, the mass 
values of the glueballs can then be converted into physical units. 

In the final stage of the calculation, glueball mass values obtained from a finite 

lattice have to be extrapolated to the continuum limit. This means that finite 
volume effects and ffnite lattice spacing errors have to be eliminated. Typi- 
cally, finite volume errors are found to be small in these calculations [6,7]. It 
is the finite lattice spacing errors that are more difficult to handle, especially 
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for the scalar glueball mass. It was observed that, the scalar glueball sector 
exhibits a dip in the extrapolation, making the extrapolation less dependable 
compared with other channels [6]. To remedy this situation, a simulation at a 
smaller values of Og has been performed. In our calculation, we also performed 
a simulation at a smaller lattice spacing where ~ 0.08/m. These two simu- 
lations now makes the extrapolation in the scalar sector more dependable and 
less sensitive to the form of the functions used in the extrapolation. 

2.1 Setting the scale 

In our simulation, the scale is set by measuring Wilson loops from which the 
static quark anti-quark potential V{R) is obtained. Using the static potential 
between quarks, we are able to determine the lattice spacing in physical units 
by measuring Tq, a pure gluonic scale determined from the static potential 
[9]. The definition of the scale ro is given by: B?dV{R)/dR\R=rf^ = 1.65. In 
physical units, ro is roughly 0.5/m which is determined by comparison with 
potential models. For a recent determination of Tq, please consult Ref. [10]. 

In order to measure the Wilson loops accurately, it is the standard procedure 
to perform single link smearing [2,6] on the spatial gauge links of the con- 
figurations. In this procedure, each spatial gauge link is replaced by a linear 
combination of the original link and its spatial staples. Each spatial gauge 
staple is weighted by a parameter Xyy relative to the original gauge link. The 
final result is then projected back into SU{3). This smearing scheme can be 
performed iteratively on the spatial links of gauge fields for as many as nw 
times. The effect of smearing is to projects out higher excitation contamina- 
tions from the Wilson loop measurements. Then, Wilson loops arc constructed 
using these smeared links. For a Wilson loop of size RxT, it is fitted against: 

M/(i?,T)^~°°Z(i?)e-^(^)^ . (2) 

Therefore, by extracting the effective mass plateau at large temporal sepa- 
ration, the static quark potential V{R) is obtained. We have measured the 
Wilson loops along different lattice axis. It is seen that the measured data 
points for the quark anti-quark potential along different lattice axis lie on a 
universal line which is an indication that the improved action restores the ro- 
tational symmetry quite well. The smearing parameter Xw used in this process 
are also hsted in Table. 1. Using smearing, we have been able to obtain descent 
plateau in the effective mass and the potential atV{R) is thus determined for 
various values of R/as- 

The static quark potential is fitted according to a Coulomb term plus a linear 
confining potential which is known to work well at these lattice spacings [7]. 
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The potential is parameterized as: 

V{R) = Vo + ejR + aR 



(3) 



From this and the definition of ro, it follows that: 



/1.65 + ec , - 

ro/as = \ 5 — • (4) 



To convert the measured result to ro/ag, we have also used the value of ^ 
taken as the bare value. As explained earlier, the renormalization effects for 
this parameter is small. The result of the spatial lattice spacing in physical 
units is also included in Table. 1. The errors for the ratio ro/cts are obtained 
by blocking the whole data set into smaller blocks and obtain the error from 
different blocks. 



2.2 Gluehall mass measurement 



To obtain glueball mass values, it is necessary to constructed glueball opera- 
tors in various symmetry sectors of interests. On a lattice, the full rotational 
symmetry is broken down to only cubic symmetry, which is a finite point 
group with 5 irreducible representations They are labeled as : Ai, A2, E, 
Ti and T2. The first two irreducible representations are one- dimensional. The 
third is two-dimensional while the last two are three-dimensional. In practice, 
one is interested in the scalar, tensor and vector glueballs in the continuum 
limit. The correspondence is the following: scalar glueball is in the Ai repre- 
sentation of the cubic group; tensor glueball is in representation E + T2, which 
forms a 5-dimensional representation; vector glueballs are in the representa- 
tion Ti. Of course, this correspondence is not one to one but infinite to one. 
Therefore, what we can measure is the lowest excitation in the corresponding 
representation of the cubic group [11,3]. 

Glueball correlation functions are notoriously difficult to measure. They die 
out so quickly and it is very difficult to get a clear signal. In order to enhance 
the signal of glueball correlation functions, smearing and fuzzying have to 
be performed on spatial links of the gauge fields [2,6,7]. These techniques 
greatly enhance the overlap of the glueball states and thus provide possibility 
of measuring the mass values. In our simulations, after performing single link 
smearing and double link fuzzying on spatial links, we first construct raw 

^ Since we only consider the glueball states with zero momentum, it suffices to 
study the cubic group. For glueball states with non-vanishing momenta, the whole 
space group has to be considered. 
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Fig. 1. The Wilson loop shapes used in constructing the glueball operators. 

operators at a given time slice t\ \vit^ = Ex^x,t}) where W^^l are closed 
Wilson loops originates from a given lattice point x = (t, x). The loop-shapes 
studied in this calculation are shown in Fig.l. Then, elements from the cubic 
group is applied to these raw operators and the resulting set of loops now forms 
a basis for a representation of the cubic group. Suitable linear combinations 
of these operators are constructed to form a basis for a particular irreducible 
representation of interest [3]. We denote these operators as {0^^{t)} where it! 
labels a specific irreducible representation and a labels different operators at 
a given time slice t. In order to maximize the overlap with one glueball state, 
we construct a glueball operator Q^^\t) = Ea "^i^^^i'^H^)) where 0^\t) = 
0^a\t) - (0 I 0^^\t) I ). The coefficients 

yCR) Qj-Q iq ^QQ determined from a 
variational calculation. To do this, we construct the correlation matrix: 

Ca/3W=E(0|(5i^n^ + T)Oi^nr)|0) . (5) 

T 

The coefficients v^'' are chosen such that they minimize the effective mass 
^effitc) = -— log 



{R)JR) 



(6) 



where tc is time separation for the optimization. In our simulation tc — 1 is 
taken. If we denote the optimal values of f^^^ by a column vector v'^''^^ this 
minimization is equivalent to the following eigenvalue problem: 

C{tc)-v^^^ ^e-'^"'^ff^'^^C{0)-v^^^ . (7) 

( R) 

The eigenvector Vq with the lowest effective mass then yields the coefficients 
Vq^^ for the operator QQ^\t) which best overlaps the lowest lying glueball in 
the channel with symmetry R. Higher-mass eigenvectors of this equation will 
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then overlap predominantly with excited glueball states of a given symmetry 
channel. 

With these techniques, the glueball mass values are obtained in lattice units 
Table 2 

Glueball mass estimates for the symmetry channel Af~^, E~^^ and at various 
lattice spacings. The entries corresponding to the highest beta value are the val- 
ues after the infinite volume extrapolation. The last row tabulated the continuum 
extrapolated result of the glueball mass values in units of l/ro- 
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and the final results are listed in Table. 2. The errors are obtained by binning 
the total data sets into several blocks and doing jackknife on the blocks. 

2.3 Extrapolation to the continuum limit 

As has been mentioned, finite volume errors are eliminated by performing 
simulations at the same lattice spacing but different physical volumes. This 
also helps to purge away the possible toleron states whose energy are sensitive 
to the size of the volume. A simulation at a larger volume is done for the 
smallest lattice spacing in our calculation. We found that the mass of the 
scalar glueball remains unchanged when the size of the volume is increased. 
The mass of the tensor glueball seems to be affected, which is consistent with 
the known result that tensor glueballs have a rather large size and therefore 
feel the finiteness of the volume more heavily. The infinite volume is obtained 
by extrapolating the finite volume results using the relation [12]: 

atM(^)(L,) = a^M^^) (oo) (l - A^^) exp(-V3^/2)/z) , (8) 

where z — Using the results for the mass of the and T2 

glueballs on 8^24 and 10^30 lattices for the same value of /3, the final result for 
the mass of these glueball states are obtained. Glueball mass values for other 
symmetry sectors are not so sensitive to the finite volume effects. Therefore, 
in Table. 2, only the extrapolated values for the smallest physical volume are 
tabulated. Other entries are obtained from 8^24 lattice results. 
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As for the finite lattice spacing errors, special attention is paid to the scalar 
glueball sector where the continuum limit extrapolation was known to have 
problems. Due to the simulation points at small lattice spacings, around 0.1 fm 
and below, the ambiguity in this extrapolation is greatly reduced. We have 
tried to extrapolate the result using different formula suggested in Ref . [7] , the 
extrapolated results are all consistent within statistical errors. For definiteness, 
we take the simple form: 

roMcias) = roMdO) + c^{as/rof + C2(a,/ro)^ , (9) 
and the result is illustrated in Fig. 2. The final extrapolated results for the glue- 
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Fig. 2. The continuum limit extrapolation of glueball mass values in scalar and 
tensor channels. The solid symbols are results from this calculation with the cor- 
responding continuum limit extrapolation represented by the solid lines. For com- 
parison, the corresponding results from [7] are also shown with open symbols and 
dashed lines. 

ball mass values are also hsted in Table. 2. The data points from our simulation 
results are shown with solid symbols and the corresponding extrapolations are 
plotted as solid lines. It is also noticed that the extrapolated mass values for 
E^'^ and T^^ channels coincide within statistical errors, indicating that in 
the continuum limit, they form the tensor representation of the rotational 
group. For comparison, simulation results from [7] are also shown with open 
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symbols and the corresponding extrapolation are represented by the dashed 
lines. We also tried to extrapolate linearly in {cIs/tqY using three data points 
with smallest lattice spacing. The results are statistically consistent with the 
results using the extrapolation (9) within errors. It is seen that, due to data 
points at lattice spacings around 0.1/m and below, the uncertainties in the 
extrapolation for the glueball mass values are reduced. 

Finally, to convert our simulation results on glueball masses into physical 
units, we use the result Tq^ — AlOMeV. The errors for the hadronic scale 
ro is neglected. For the scalar glueball we obtain Mg(0"'"+) = 1730(90)MeV^. 
For the tensor glueball mass in the continuum, we combine the results for the 
T2++ and E++ channels and obtain Mg(2++) = 2400(95)Mel^ for the tensor 
glueball mass. 



3 Discussions and Conclusions 

We have studied the glueball spectrum at zero momentum in the pure SU (3) 
gauge theory using Monte Carlo simulations on asymmetric lattices with the 
lattice spacing in the spatial directions ranging from 0.08/m to 0.25 fm. This 
helps to make extrapolations to the continuum limit with more confidence for 
the scalar and tensor glueball states. The mass values of the glueballs are con- 
verted to physical units in terms of the hadronic scale Tq. We obtain the mass 
for the scalar glueball and tensor glueball to be: 771^(0++) = 1730(90)Mey 
and 777.0(2++) = 2400(95)Mey. It is interesting to note that, around these 
two mass values, experimental glueball candidates exist. Of course, in order 
to compare with the experiments other issues like the mixing effects and the 
the effects of quenching have to be studied. 
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